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The tacnode Riemann-Hilbert problem 

Arno Kuijlaars* 



Abstract 



■^r , The tacnode Riemann-Hilbert problem is a 4 x 4 matrix valued 

r ) ' RH problem that appears in the description of the local behavior of 

two touching groups of non- intersecting Brownian motions. The same 
RH problem was also found by Duits and Geudens to describe a new 
2 . critical regime in the two-matrix model. 

Delvaux gave integral representations for some of the entries of 
the 4x4 matrix. We complement this work by presenting integral 
representations for all of the entries. As a consequence we give an 
^ ', explicit formula for the Duits-Geudens critical kernel 

l> 

CS . 1 Introduction 

'nT . The tacnode Riemann-Hilbert (RH) problem is a 4 x 4 matrix-valued RH 

ir I problem that first appeared in the asymptotic analysis of two touching 

groups of nonintersecting Brownian motions, a so-called tacnode. The po- 
sitions of the non-intersecting Brownian motions are a determinantal point 
process that in a double scaling limit around the tacnode leads to the tac- 
node process. The one-time correlation functions of the tacnode process 
^_' were expressed in terms of the tacnode RH problem in [T2]. The tacnode 

RH problem is related to the Hastings-McLeod solution of the Painleve H 
equation as was also discussed in [12]. 

The tacnode problem was also analyzed in [21 [161 122] using different 
techniques. In these papers the tacnode kernel and its multi-time extension 
are expressed in terms of integrals with resolvents of Airy integral operators 
acting on a half-line. These expressions are very different from the RH 
formulation. 
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Recently, Delvaux |10j made the connection between the two sets of 
formulas by presenting integral representations for some of the entries of 
the solution of the tacnode RH problem. These entries are exactly the ones 
that play a role for the tacnode kernel in |12j. With these explicit formulas 
Delvaux could make the connection between the formulas in [12j and the ones 
by Ferrari and Veto [16] for the asymmetric tacnode. The paper [lOj was 
inspired by the paper |6j by Baik, Liechty, and Schehr, where a connection 
between different sets of formulas for the maximal height and position of 
the Airy2 process was made. 

The aim of this paper is to complement the work of [10] by providing 
integral representations for all the entries of the tacnode RH problem. Some 
of these remaining entries appear in the description of a critical kernel ap- 
pearing in the two-matrix model as shown by Duits and Geudens [13]. We 
therefore find explicit integral formulas for the Duits-Geudens critical kernel. 

In section [2] we recall the tacnode RH problem with some of its prop- 
erties, and in particular the connection with the Hastings-McLeod solution 
of Painleve H. The main results of this paper are stated in Theorems 12.51 
and 12.81 below. We compare the solution of the tacnode RH problem with 
the explicit solution of the usual 2x2 matrix-valued RH problem for the 
Hastings-McLeod solution in section 12.61 

The proofs of the results are in section [3l A key role is played by Lemma 
13.21 that describes solutions to a certain ODE system (13. 2p . The proof of 
this lemma follows along the lines of certain proofs in [10]. We give full 
details about the calculations in section O Following [10] we briefly mention 
the tacnode kernel in section 14.11 The implications of Theorem 12.81 for the 
Duits-Geudens critical kernel are discussed in section [H 

A variation of the tacnode RH problem for the hard-edge tacnode and 
the chiral two-matrix model appears in [9l [11]. It may be possible that 
explicit integral representations for the solution of these RH problems can 
be found as well. Other recent contributions [H O [71 [19] discuss further 
connections and properties of the tacnode process. 

2 Statement of results 

2.1 The tacnode RH problem 

The tacnode RH problem asks for a 4 x 4 matrix-valued function 

M : C \ Fa/ ^ C^""^ 
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Figure 1: The figure shows the jump contours F^ in the complex plane 
and the corresponding jump matrices J^, A: = 0, ... ,5, in the tacnode RH 
problem. Qk is the sector bounded by the two rays r^ and r^+i. 



which is defined and analytic outside a set Tm which is a union of six rays 



M 



k=0 



Tk = {z\ argz 



|.). 



(2.1) 



as shown in Figure [TJ Each ray is oriented from the origin to infinity. The 
orientation induces it-sides on each ray, where the +-side is on the left and 
the — side is on the right as one traverses the ray according to its orientation. 
We ask that M has continuous boundary values M± on each of the rays that 
satisfy the jump condition 



M. 



M-Jk on Ffc for fc = 0, 1, 



(2.2) 



where the jump matrix J^ on F^ is also shown in Figure [H 

The RH problem depends on a number of parameters ^i, r2, si, S2, t that 



appear in the asymptotic condition for M via two functions 

Oi{z) = ln{-zf^ + 2si{-zf/\ z G C \ [0, oo), 

(2.3) 
e2{z) = \r2z'^'^ + 2s2^^/^ zeC\ (-00, 0], 

where we use principal branches for the fractional exponents. The asymp- 
totic condition is 

M{z) ={l + M^ + 0{z-')\ diag((-^)-l/^ z-^'\ {-zf'\ z'") 



1 

^ —F 
y/2 



/I -i 0\ 

1 i 

-i 1 

V i l/ 



diag ( e-^l(^)+^^ e-'2(z)-rz^^e,{z)+rz^ ^B,{z)-rz 



(2.4) 
The residue matrix M^^> is independent of z but depends on the parameters. 



Remark 2.1. In the papers [121113] the tacnode RH problem was formulated 
on a union of ten rays. Here we choose to combine the two jumps in each 
of the open quadrants which reduces the number of rays by four. It is easy 
to see that the two RH problems are equivalent. 

Proposition 2.2. Suppose that the parameters ri,r2, si, S2,t are real with 
ri > and r2 > 0. Then the RH problem (J2.2p - (j2.4p has a unique solution. 



Proposition 12.21 was proved in the case r = by Delvaux, Kuijlaars, and 
Zhang [12j, and in the case ri = r2 = 1, si = S2 with general r by Duits 
and Geudens [13] . The proof in [13] extends to the general case as noted by 
Delvaux [10]. 



2.2 The Painleve II connection 

The tacnode RH problem is related to the Hastings-McLeod solution of the 
Painleve II equation, as was noted in the cited papers [ini IS [13]. The 
Painleve II equation is 

g" = tq + 2g3 (2.5) 



and the Hastings-McLeod solution is the unique solution of (|2.5p that sat- 
isfies 

q{t) = k\{t){l + o{l)) ast-^+00 



where Ai denotes the usual Airy function. We also need the related function 

u = ci'^ -tq^ -q^, (2.6) 

which satisfies u' = —q^. 

Proposition 2.3. The solution M of the tacnode RH problem satisfies a 
differential equation 

—M = UM, (2.7) 

oz 

with a matrix U that is explicitly given in terms of the Hasting s-McLeod 

solution q of Painleve II and u from, (|2.6p as follows: 

U = 



iri 








ir2 


,2 " 


Vni 



(2. 



iJere ri,r2, si, S2)''" a'^e the parameters in the problem, 

C=irf + r,^)^/^ (2.9) 

T^-g^^'-^^^^^). (-0) 

and the Painleve functions q,q' , and u that appear in ()2.8p are evaluated in 

2 f si S2 2r2 \ , ^ 

t = - - + --^ 2 • 2.11 

C Vn r2 rf +r|y 

Proof. See Delvaux ^lOj. Note however that the notation in ^lOj is slightly 
different from ours. The constant r used in [1^ is equal to ^^ ^ r, the 



T-i 



constant D in [TO] is equal to ^7^7 and a is used instead of t. D 



It is known that the Hastings-McLeod solution of Painleve II has no 
poles on the real line [2D]. Therefore the linear system (|2.8|) is well-defined 
for every choice of real parameters ri,r2,si,S2, and r. 

There are also linear differential equations 

r\ O O 

-^—M = ViM, ^—M = V2M, —M = WM, (2.12) 

OSi OS2 OT 



with explicitly known matrices Vi, V2 and W. The Painleve II equation 
(|2.5p arises as the compatibility condition for ()2.7p and (j2.12p . which can be 
viewed as a Lax pair. Note that the usual Lax pair for Painleve II is of size 
2x2, see [171 [1^, ^^^ also section \TQ\ below. 

From (|2.7p it follows that each column of M is a solution to the linear 
system of ODEs 

fjTfy 'T' 

— — = Um, m = (mi 1712 Tn^ 1714) (2.13) 

oz ' 

with V given by (12. Sp . To specify a solution of (J2.13P it is enough to give 
m\ and r?i2 since 



irims = mi - [T - Si + -^jmi ■=p::;m2, (2-14) 



, \/fiQ 2 u ^ 

ir2m4 = ^2 + 7 rn-i - (-r + S2 - -^)m2, (2.15) 

-v/r2C C 



which follows from the special structure of U in (j2.8p . 

2.3 Tracy- Widom functions 

We are going to construct six solutions m,^^', j = 0, 1, . . . , 5, of (|2.13p . The 
formulas are based on functions that were first introduced by Tracy and 
Widom [MIES]. 

For t G R, we use Kt to denote the integral operator on [t, 00) with the 
Airy kernel 

Ai(x)Ai'(y)-Ai'(2;)Ai(y) 



x-y 
Thus Kt : L'^{[t, 00)) -^ L'^{[t, 00)) is defined by 



(2.16) 



W)ix)=l mAM^AM^Mf^y^,y. (2.17) 



It is known that I—Kt is invertible on L/^ilt, 00)), and we define two functions 
by 



We also put 

Rt{x,y) 



Qt = {I- Kt)-' Ai, 


(2.18) 


Pt = il- Kt)-' Ai' . 


(2.19) 


Qt{x)Pt{y) - Pt{x)Qt{y) 


(2.20) 


x-y 
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Both Qt and Pt are continuous (in fact real analytic) functions on [t, oo). 
It is known that 

Qt{t)=q{t), Rt{t,t)=u{t), (2.21) 

where q is the Hastings-McLeod solution of (|2.5|) and u is given by (|2.6|) . 
see e.g. section 2.3 in [25j and in particular formula (25). 

Lemma 2.4. Both Qt and Pt extend to entire functions on the complex 
plane and 

Q,(x) = (1 + 0{x-^/^)) Ai(x) = O (x-i/^e-i^''") , (2.22) 

p^{x) = (1 + 0{x-^)) Ai'(x) = O (x^/^e-i^''") , (2.23) 

as X —7- oo, uniformly for — vr + e < argx < tt — e, for every e > 0. 

Proof. According to the Its-Izergin-Korepin-Slavnov theory on integrable 
operators ^1] (see also [8]) we have that 

where y is the unique solution of the RH problem: 

• y : C \ [t, oo) ^ C^^^ is analytic, 

• Y has continuous boundary values y+ and Y- on (t, oo) that satisfy 

y+(x) = y_(x)(/-2^^(^^^^^)^)(Ai'(x) -Ai(x))), (2.25) 

for X G (i,oo), 

• y(x) = 0(log |x — t|) as X — )■ t, 

• y(x) = / + 0(l/x) as X — )• oo. 

Thus X G C I— 7- y(x) I . ., , \ I provides the extension of ( _ | into an 
\ki {x)J yPtj 

entire function on the complex plane. The lemma then follows because of 

the asymptotic behavior of Y and the well known behavior 

AKx) = ^-^e-i^^^^(l + 0(x-/^)) 



(2.26) 
as X — )• oo, I arg x\ < vr — e of the Airy function. D 



Ai'(x) = ^^e-3- ' (l + 0(x-3/2)) 



By ()2.20p we then also have that Rt extends to an entire function in the 
complex plane and for every fixed y, 

i?t(x,y) = 0(x-3/V^''') (2.27) 

as X — 7- oo with — vr + e < arg x < vr — e for some e > 0. 

2.4 Six solutions of f l27[3|) 

The functions Qt{x) and Rt{x,t) from ()2.18p and ()2.20p appear explicitly in 
the integral formulas we have for the solutions of ()2.13p . We use 

We also put 

^ = ^-^^^ /^ = ;:^-72^ (2-28) 

'l "T '2 'l "T '2 

and we recall the definitions of C, 7 and t in ([23]), (l2lil and (IZTTT) . 

Theorem 2.5. There are six solutions m^^', j = 0, . . . ,5 of (j2.13p whose 
first and second components are given below. In all cases the third and fourth 
components are as in (J2.14p - (j2.15p . 

The six solutions are given as follows. 

Solution m(°) Let Fo{z) = Ai (r^^^z + ^j e-'^^^'' . Then 

(mf\z)\ ^ ^^i/e^x. ( -0t r Fo{z + C{x - t))Q,{x)dx \ 
\mf\z)) ' \F^{z) + j^ Fo{z + C{x - t))Rt{x,t)dx) ' 

(2.29) 

Solution m(^) Let Fi{z) =ujAi(uj ( r^z + ^ ) ) e'^'^f'^ Then 

(m['\z)\ _ /- 1/6 xz fF^i-z)+JT"F,i-z + C{x-t))Rt{x,t)d: 



l-Kr-i 

(2.30) 



■x 



A'\^) ^1^1 _^^J- F^(_, + c(x-t))Q,(x)dx 



Solution m(2) Let F2{z) = u? Ai iu? (r^^^z + ^ J J e'''^^' . Then 



,m^2\z)J ' \F2{z) + j;;^'^F2{z + Cix-t))Rt{x,t)dx^ 

(2.31) 



Solution m(3) Let F^{z) = Ai (rl'^z + ^ j e'^'^^^ Then 

(mt\z)\ _ r- 1/6 A. /i^3(-^) + r^3(-^ + ^(x-t))i?t(:E,t)dx\ 

(2.32) 

Solution mW Let F^[z) = uj Kiiu ir^J^z + ^ ) ) £"'^2'^^ T/ien 

K^('^U^//V^f -^/r^F4(z + C7(x-t))Q,(x)dx \ 

W^(^)y ' v^4(^)+/r"'^4(^+c(x-t))i?i(x,t)dxy 

(2.33) 
Solution m(5) Lei ^5(2) = w^ Ai iuJ^ (rl^^z + ^ ) ) 6"^'?''^ Then 

(^?\A _ nr 1/6 A. /i^5(-^) + r"^5(-^ + C^(x-t))i?i(x,t)dx^ 
U^)(z)J-^'"'^ ' 1^ -l^jrF^i-^ + C{x-t))Q,[x)dx ^ 

(2.34) 

The proof of Theorem 12.51 is given in section 13.31 

Remark 2.6. The solutions m^' and tti'^^ were found by Delvaux [10] , 
These are the solutions for which the integrals are taken over the real interval 
(t, cxd). The other solutions are new and their identification is the main result 
of this paper. 

Remark 2.7. The integrals in ^(T^ and (fOil) start at t G M and end 
at infinity at asymptotic angle 27r/3. Similary, the integrals in (|2.30p and 
(|2.33p end at asymptotic angle — 27r/3. 



By Lemma 12.41 we have that Qt{x) and Rt{x,t) are entire functions in 
X. To see that the integrals in (j2.3ip indeed converge, we note that both 
Qt{x) = 0(e-i^'^') and Rt{x,t) = 0{e-h^^') by (12:22]) and (E^ZD, and so 

as a; — 7- oo with argx = 27r/3. Since by ()2.9p 

the integrands decay at an exponential rate at infinity, and the integrals in 
(|2.3ip converge. 

Similarly, the integrals in (j2.30p . (j2.33p and (|2.34p converge. 

2.5 The solution of the tacnode RH problem 

The vector m^^', j = 0, . . . , 5 turns out to be the recessive solution of (I2.13P 
in the sector 

Sj={z£C\ -f +if < argz < f + jf}. (2.35) 

Note that Sj is the sector of angular width tt/3 with Tj as its bisector, see 
p.ip . The constant prefactors in the definitions (I2.29p -( r2.34p are chosen 
such that m^^> appears as one of the columns of M in the two sectors ^j-i 
and Qj that intersect Sj, where it is understood that ri_i = ^5. As such 
the vectors m^^' are the building blocks for the solution of the tacnode RH 
problem. 

Our main result is the following. 

Theorem 2.8. The solution of the tacnode RH problem (|2.2p - (j2.4p is given 
by 

M = (m(3) m(°) m(i) m^^)) 

in the sector tt/3 < arg z < 2tt/3 around the positive imaginary axis. In the 
other sectors it can be found by following the jumps (12. 2p in the RH problem, 
and by using the non-trivial relations 

^(0) ^ ^(3) ^ ^(1) _ ^(5) ^ ^(2) _ ^(4) ^2.36) 

among the six solutions m*-°\ . . . , m*-^) of (j2.13p . 

Explicit expressions for M in all sectors ilj are given in Figure \M 

The proof of Theorem 12.81 is given in section 13.41 

The relations (I2.36P among the recessive solutions of (12.13P are quite 
remarkable as they do not follow in a straightforward way from the integral 
representations. It is an open problem how to prove these identities in a 
direct way from the formulas (I2.29p - (l2.34p . 
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(m(3) mW m(i) m^^)) 
il2 \ Q / ilo 




o\ 



to(i) m(0) m(5) m^-*) - m^^)! 






Figure 2: Structure of the solution of the tacnode RH problem. The column 
vector 771"-* is the recessive solution of (I2.13P in the sector Sj. 



2.6 Comparison with the 2x2 RH problem for Painleve II 

The usual RH problem for Painleve II is of size 2x2 with a contour Tq, as 
in Figure El Then ^ : C \ F^ — ^ C^^^ is analytic satisfying *+ = ^_ J^ on 
F^ with jump matrices J^ that are also shown in Figure lU and with the 
asymptotic condition 

^(z) = (/ + 0(l/z))e-^(t^'+*')"3 , as z ^ 00 

where (^3 = \ p. -, I • The Stokes multipliers si, S2, S3 that appear in the 

jump matrices are assumed to satisfy S1S2S3 + si + S2 + S3 = 0. Then a 
unique solution to the RH problem exists (except for an at most countable 
number of t values), and $ satisfies the differential equations (Lax pair) 

d (-Aiz^ - i{t + 2g2) Azq + 2iq' \ 

5; V ^zq-2iq' Aiz^+i{t + 2q^))^ ^^""^^^ 

d (-iz q\ ^ 

at \ q izj 

where q = q{t) is a solution of the Painleve II equation, determined by the 
Stokes multipliers, see e.g. pTlfTS]. 
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Figure 3: Jump matrices in the 2x2 RH problem for Painleve II. The Stokes 
multiphers satisfy S1S2S3 + si + S2 + S3 = 0. The Hastings McLeod solution 
corresponds to si = 1, S2 = and S3 = —1. 



The special case si = 1, S2 = 0, S3 = —1 leads to the Hastings McLeod 
solution of the Painleve II equation. In this case there is an explicit formula 
for the solution of the RH problem, which is contained in the recent work 
of Baik, Liechty, and Schehr [6j , see also [5j . The solution is built out of the 
two column vectors 






4 3+,,)/l + /^°°e-2^(--*)^i?i(x,t)dx 



-^(. 



^i{-^z'-^+tz) 



It 



00 -2i(x-t)z 



Qt{x)dx 



l + J^'^e^'^^-'>Rt{x,t)dx 



(2.38) 
(2.39) 



that satisfy the vector analogue of ()2.37p . Then the solution of the RH 
problem is 

^(z) = {il^W il^i^)) 



in the two sectors — ^ 



< arg z < ^ and ^ < arg z < ^. The solution in all 



sectors is given in Figure [H 
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; 1) (^^^^ ^^^'+^^^^) (; "i' 



Figure 4: Solution of the 2x2 RH problem for the Hastings-McLeod solution 
of Painleve II in terms of the columnn vectors V' and V' • 



Comparing this solution with the explicit solution M of the 4x4 matrix- 
valued RH problem, we see that the same functions Qt{x) and Rt{x,t) ap- 
pear in the solutions, but there does not seem to be a direct way to go from 
^ to M. The 4x4 RH problem and its associated Lax pair therefore pro- 
vide a genuinely different characterization of the Hastings-McLeod solution 
of Painleve II. 

3 Proofs of Theorems [23] and [2781 

3.1 Transformation to second order system 

It will be convenient to transform the first order system (j2.13p to a second 
order system. The transformation also removes the parameter 7. 

Lemma 3.1. Suppose m= (mi 1112 m^ m.4) satisfies (I2.13p . Then 



V'(-) = n'!'l)=«"'M''^"'^'M (3-1) 
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satisfies 













-4' 

dz 


+ SV' 




(3.2) 


A 


( 2rl^, 


1 -'. 












(3.3) 


B 


f-rlz + 2ris 


i + Cr 
-Crh' 


iV- 


4 2 
r^fj, 

2 

^2 


■z + 2r2S2 + Cr'^q'^ - 


„4,,2 


) • (3.4) 




Conversely, 


ifip 


= {i'i, 


V'sf 


satisfies 


(13.21) with A and B 


given 


6y(I23D- 


13.4 


1 then m = 


(mi, 


m2, ma, 7714)' 


^ w/iere 









mi(z) = e^^7"'^V^i(^), "12(2) = e^'M^) (3.5) 



and ma, m4 are given by (I2.14p - (12.15p satisfies (12.13p . 

Proof. This is a straightforward calculation and we will not give full details, 
see also [TUl Proposition 2.12]. 

Let us just note that for general A one obtains ()3.2p with 

/2r - 2A CVfg 
^ - ' -CV^g -2t - 2A 



-rjz + 2risi + Crfq'^ - (r - A)^ -Crfg' + ((rf - r2)T + (rf + r2)A))^' 

-Grig' - ((rf - ri)T + (r^ + r^A))^ r^z + 2r2.2 + Cr^q^ - (r + A)^ '' ^ 

Because of the choice (|2.28p for A we have that the off-diagonal entries of B 
simplify since the terms with q disappear. By (j2.28p it is easy to check that 

r — A = rffi, r + A = r|/i 

and we obtain iKM and (1331). □ 




3.2 Solutions to the second order system (13.21) 

We denote 00 = e^'^*'^ as before, and we let 

yo{x) = Ai(x), 2/1(2;) = ujAi{ujx), y2{x) = uJ^ A\{uj'^x) 

be three solutions of the Airy differential equation y" = xy. 

The following lemma gives solutions to the second order system (13. 2p . 

14 



Lemma 3.2. For fc = 0, 1, 2, we put 



F{z) = Vk 1 rll'^z + -jj^ J exp(-r2Ai2;). (3.6) 



Then the vector li" = ( , 1 where 



Mz) = - F{z + C{x-t))Qt{x)dx, (3.7) 

F{z + C{x-t))Rt{x,t)dx, (3.8) 

is a solution of ()3.2p . 

Proof. For A; = 0, this is proved in |10t section 5.2], although with somewhat 
different notation. The proof uses the differential equation for F 

F"{z) = -2rliiF'{z) + {rlz + 2r2S2 - 4fi^)Fiz), (3.9) 

see also [101 Equ. (5.2)] (where F is called bz{x) and /i is called r), and the 
fact that the integrals in (j3.7p and (j3.8p converge for A; = 0. We already 
noted, see Remark 1.6, that the integrals in (13. 7p and (13. 8j) converge for 
k = 1,2 as well. 

What is also used in [10] are a number of identities for the functions Qt, 
Pt and Rt introduced in ()2.18p - (l2.20p . namely the differential identities (see 
also iMlES] or [i section 3.8]) 

^^Rt{x,y) = -Rt{x,t)Rt{t,y), (3.10) 

^^Qt{x) = -Rt{x,t)Qt{x), (3.11) 

^^Pt{x) = -Rt{x,t)Pt{x), (3.12) 

and identities for the x-derivatives of Qt and Pj, 

Q[{x) = Pt{x) + q{t)Rt{x, t) - u{t)Qt{x), (3.13) 

P;(x) = xQtix) + p{t)Rt{x, t) + u{t)Pt{x) - 2v{t)Qt{x). (3.14) 

The identities ()3.10p ~ p.l4p of course extend into the complex plane. 
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Also for k = 1,2, we have that F satisfies (|3.9p . Since the integrals 
in dSZD and ([33]) converge for A; = 1,2, and the identities (f3l^ - ([3J2D 
and (|3.13p - (|3.14p remain vahd for x in the complex plane, we can follow the 
proof in |10l section 5.2], making proper modifications due to some change in 
notation. For convenience of the reader we provide the detailed calculations 
in the appendix (section [5]) using the notations of this paper. D 

3.3 Proof of Theorem [231 
We can now prove Theorem 12.51 



Proof. In view of Lemmas 13.1 1 and l3. 21 we find three solutions of ()2.13p . After 
multiplication by appropriate constants these are the solutions m^^', m^"^' 
and m('') given by the formulas (12:29]) . (lOTT) . (12:33]) in Theorem [231 

The other solutions follow from a symmetry in the system (|2.13p . Namely, 
if m{z) is a solution of (I2.13P then 

J 0\ , . , /O 1 

m[—z), J 



,0 -J) ^ " V Oy 

solves (I2.13P as well, but with the change of parameters 

ri o r2, si ^ S2, T I-?- T. (3.15) 

The constants C, t and n do not change under this change of parameters, 
but 

7 i-^ 1~^: A H> -A, 

see the formulas ([2:9]) - ([2lT]) and ([228]). Thus one solution of (|2T3]) leads 
to another by a change of sign z i— )• — z, a change of parameters ()3.15p . 
combined with an interchange mi o m2- In this way the solutions m^^', 
mP'' and m^'^> lead to the solutions m^^\ m^^', and m^^'. 

This completes the proof of Theorem 12.51 D 

3.4 Proof of Theorem [278] 

Proof. The column vectors of M are in each sector ^^ a basis of the vec- 
tor space of solutions of (|2.13p . In each sector the four different columns 
represent the four different types of asymptotic behavior, as given by ()2.4p . 
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Namely, if e^ denotes the kth unit vector, 

M(z)ei ~ e-^i(^)+^^ 



(3.16) 

M(z)e4 ~ e' 



92(2:) — T2 



as z ^ 00. 

Distinguished solutions of (|2.13p are those solutions that are recessive 
in a certain sector, i.e., they are smallest possible as z — )■ 00 in that sector. 
Recessive solutions are unique up to a multiplicative constant. From (|3.16|) 
and ()2.3p one sees that M{z)ei is a recessive solution of ()2.13p in the sector 
S3, M{z)e2 is a recessive solution in 5*0, M{z)e3 is a recessive solution in 
sectors Si and ^5 and M{z)e4, is recessive in sectors ^2 and 5*4. 

It will turn out that m^^> is the recessive solution of (12.13P in sector Sj. 
This then implies, for example, that m^^' is a multiple of M[z)e2 in Sq. 
The constant \/27rr2 in the definition ()2.29p of m^^' has been chosen so 
that M{z)e2 = rn^^' in 5*0. Since Sq has a non-empty intersection with both 
Qq and r^s it then follows that m^^' appears as the second column of M in 
sectors Qq and Q5. 

Similarly, m'^^' is in the third column of M in the sectors Qq and fii, 
m^'^' is in the fourth column of M in Ui and ^2, and so on. 

This leads to the partial solution of the RH problem given in Figure [5j 

Then we complete the solution by following the effect of the jump ma- 
trices Jfe, see Figure [H and this leads to the full solution of the tacnode RH 
problem as given in Figure [2j 

We prove in more detail that m^^' is the recessive solution of (I2.13P in 
the sector S^^' and that m^^) = Me2 in rig U ^5. The arguments for the 
other solutions m^^> can be done in a similar way. 

The asymptotic condition (|2.4p tells us that 



as z — )• 00. It will be enough to show that 7712 has the same asymptotic 
behavior as z — )• 00 in the sector ^o in order to conclude that m-^' = Me2 
in ilo U ^5. 
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im^"' * * 771 




[m^ ' * 



* to'°^ m*^^' *J 



Figure 5: Partial solution of the tacnode RH problem with only those 
columns that are recessive solutions of ()2.13p in parts of certain sectors. 
The * entries denote columns that still have to be determined. 



Recall that Trig is given in ()2.29p . Since A — r^/i = — r we can write 



,2/3^ , 2S2 



+ V2^rl^^e-^' r Ai (rl'^{z + C{x - t)) + ^\ e-''2^^(^-*)i?,(x, t)dx. 

(3.18) 



From the asymptotic behavior of the Airy function, see (I2.26p . we have as 
z — )• oo with I arg z\ < vr — e, 



(see (I23D for 62), and 

Aifr2'/'(z + C7(x-t)) + ^ 
V '■2 

with a 0(2^"*^'^) term that is uniform for x > t in case | argz| < 7r/6. Thus 



2S2_ 
1/3 
'2 



1 ^_e,(.)-r,C(x-t).V2 / ^ ^ _i/2. 

2V^r2'/'zi/4 V 
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the first term in the right-hand side of ()3.18p is 

as z — )• oo, while the second term is (for z — )• oo in So) 

f-OO 



As z — )• oo in So, the main contribution to the integral in (|3.21|) comes from 
the endpoint x = t. K crude form of Laplace's method, see e.g. [23], shows 
that the integral is 0{z~^''^) and therefore we find that the second term in 
the right-hand side of (fXTSD is 0(z-3/4)g-e2(2)-Tz_ Together with ^^?M 
this gives 

^(0)(^) = (^z-^l^ + 0{z-^'^)\ e-^^W— (3.22) 

as z — )• oo in S'o, which indeed agrees with (j3.17p . 

This completes the proof of Theorem 12.81 D 

Remark 3.3. A more refined application of Laplace's method would lead 
to 



V r2C y a/2 

and similarly 

where it is used that Qt{t) = Q{t), Q't{t) = pit) and Rt{t,t) = u{t). This 
leads to formulas for certain entries in M^"^' (this is the residue matrix from 
the asymptotic condition, see (|2.4p ). but we do not discuss that here. 

4 Correlation kernels 

As already mentioned in the introduction, there are two correlation kernels 
in random matrix theory that can be expressed in terms of the tacnode 
RH problem, namely the one-time correlation kernel for the tacnode process 
[12j . and a critical kernel in the two-matrix model [13]. The implications of 
the explicit form of the solution of the tacnode RH problem for the tacnode 
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process was discussed by Delvaux [lOj, who made the connection between 
[12] and the different sets of formulas derived in [21 [3l [161 122] • We briefly 
discuss it in the section 14.11 

The critical kernel in the two-matrix model is due to Duits and Geu- 
dens |13J . Theorem 12.81 yields an explicit integral representation for this 
correlation kernel, as we show in section 14.21 

4.1 Tacnode kernel 

The tacnode kernel is 



i^*"'^(x,y;ri,r2,si,S2,r) = — -^ (O 1 l)M-\y)M{x) 

lm[x — y) 



1 




(4.1) 
where M denotes the analytic continuation of the restriction of M to the 
region Oi around thepositive imaginary axis. This was obtained for r = 
in [12l Definition 2.6]j and for general r £ M in [TOl section 2.2]. Thus by 
Theorem 12.81 we have 

M = (m(3) mW m(i) m(2)) , (4.2) 

see also Figure [21 

Because of the symmetry, see [9, Lemma 5.1] or [10^ Lemma 3.1], 

o -A ,., .T f o I 



M~'iz;T)=^^ ^JM{z;-Ty ^_^ ^j (4.3) 

(we use M{z; r) to denote the dependence on r, and / is the 2x2 identity 
matrix), we can also write 

K^'"'{x,y;ri,r2,si,S2,T) 



^ '110 0)M{y;-Tf(^^ ^q) M{x;t) 



27ri{x — y) 



■ (4.4) 



\0j 



^There is a misprint in formula (2.47) in T^. It should be M ^{v)M{u) instead of 
M-'^{u)M{v). 
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Thus K^"-^ only depends on the sum of the first two columns of M. If we 
put 

(4.5) 



fh = mS '+111}' 



then we get by (fOj) and (jO]) 



1 



O I 
-I o 



fri{y\-TY' [^j t.] fh{x;T) 



1 



27ri{x — y) 



27ri{x — y) 
(-mi(x; T)m3{y, -r) + m3(x; T)mi{y, -r) 

-m2ix; T)mi{y, -r) + m^ix; T)m2{y, -r)) . (4.6) 



A remarkably simple expression for 



dx dy 



K {x, y) can be obtained 



from (|4.ip . Using the differential equation ()2.7p for M and the formula (|2.8 
for ?7, we obtain 



^ + |l^""(-^' 



1 



27rz(x — y) 



(0 1 l)M"^(y)(C/(x)-?7(y))M(x) 



/1\ 

1 





— (001 l)M-\y){nE^^i-r2Ei^2)M{x) 



1 



Vo/ 



where Ej^k is the matrix with 1 in position j, k and otherwise. Combining 
this with ()4.3p we get 



d d \ 1 

— + — j K^"-''{x, y) = Y' (''i"^i(y; -T)mi{x; r) - r2rh2{y; -r)m2(x; r)) . 



(4.7) 

where m = m^'^' + tti*^^^ as in ()4.5p . 

Delvaux [TU] further analyzed (|4.6p , (|4.7p using the formulas (|2.29p and 
(j2.32p for TTT.''^' and m^^', and showed that the expression for the tacnode 
kernel agrees with the one given by Ferrari and Veto |16] . 
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4.2 Duits-Geudens critical kernel 

The Duits-Geudens kernel appears in a critical regime in the two-matrix 
model [13], where it was obtained from a Riemann-Hilbert analysis based 
on |14t I15| . It is a remarkable fact that it can be expressed in terms of the 
solution of the tacnode RH problem (J2.2p ~ (j2.4p with the special choice of 
parameters 

ri = r2 = 1, si = S2 = s, re M. (4.8) 

With those parameters we have by ([2:91) " ([2TTT]) . and (p:28|) . 

C = 2^/'\ 7 = 1, t = 2^/3(25 -r^), A = 0, ^i = t. 
The formula for the critical kernel is 

n\ 
1 



K'^"\x,y;s,T) = ^—^-—{-l 1 O) Miix)-^M{iy) 







(4.9) 



which can easily be obtained from the formulas (2.13) and (2.15) in |13j . 
Here M is the solution of the tacnode RH problem with parameters ()4.8p . 
We can use (14.31) and Theorem 12.81 to rewrite (14.91) as 



K''"Hx,y;s,T) 

= ^ (0 1 -l)M{ix;-Tf(^^ ^^M{iy-T) J 



2'Ki{x — y) 



1 (m«-m(2))(zx,-r)^f^^ ^ Wm^ + m(3))(iy,r). 



2m{x-yy ^^ ' ^ V--^ O^ 

(4.10) 

The critical kernel thus depends on m = m^^' + m^^\ see (|4.5p . which also 
appeared in the tacnode kernel, and on 

m = m(^) - m(2) = m(^) - m^^^ (4.11) 

The second identity in (14. lip holds because of (I2.36p . 

Another formula for K^'^'^^ comes from differentiating ()4.9p with respect 
to s. There are differential equations ^j- = ViM and ^j- = V2M, which if 
s = si = S2 as in (j4.8p leads to 

= VM, V = Vi + ^2, 
as 
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and V is given explicitly in [9, Proposition 5.11]. This formula implies 



Vix) - V{y) = -2i{x - y) {E^,! + i?4,2) 

from which it follows that 
d 



(4.12) 



ds 



M-^{ix)M{iy) = M-^{ix) iV{iy) - Viix)) M{iy) 



= -2{x - y)M-\ix) (^3,1 + ^4,2) M{iy) (4.13) 

Thus from (filQ]) and (jM]) we get 
d 



ds 



K'"'{x,y;s,T) 



m 



(-1 1 0)M{ix)-\E3,i + Ei,2)M{iy) 



in 



(0 1 -l)M{ix;-Ty{Ei^i + E2,2)M{iy;T 



/1\ 
1 



('\ 

1 





(4.14) 



which is somewhat similar to the expression (j4.7p for the tacnode kernel, 
except that the two solutions (j4.5p and (j4.1ip of the ODE (j2.13p are now 
involved. Indeed, by (j4.14p and the solution of the tacnode RH problem 

Q . — 1 _ 

—K''"\x, y; s, r) = —m{ix; -t^ {En + E22)m{iy; r) 
OS 7^^ 



-1 



vrz 



[mi{ix] -T)mi{iy;T) + m2{ix; -T)m2{iy]T)] , (4.15) 



which is a rank two kernel. 

Let us check that (j4.15p is real-valued. The symmetries of the tacnode 
RH problem, see [U Lemma 5.1], 

o -M'^Ao -/ ' 



M{z) ■- 

M{-z) : 
imply for z = ix with x real 



'oZ)-Mo-j}' ^ 



1 

1 



Miix)={^ ^)Miix)(^ ^ 
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In view of the solution for M in Theorem 12.81 this means that for real x,y, 



rrii [iX) = 7712 (*^j) ''^1 (1-X) = 710,2 [ix), 



ml {iy) = 7X12 (^y)) ""^1 (%) = ^-2 iw)- 
Using this and ([43]), (|4TT]l in (jilSl) . we find 

9 —2 

—K''"\x, y; s,t) = — Im [mi(ix, s, -r)mi(iy, s, r)] (4.16) 

which is real- valued (as it should be). 

Since K'^''**(x, y; s, r) — ;■ as s — )■ +00, we recover i^'="* from (I4.16P after 
integration with respect to s 

2 f°° 

K''"\x,y;s,T) = - Im [mi(ix, s', -r)mi(iy, s',r)] ds' (4.17) 

which is maybe the simplest form for the Duits-Geudens critical kernel. 

5 Appendix: Proof of Lemma 13.21 

Proof. Throughout the proof of Lemma 13.21 we simply write 00 instead of 
oow'^'^, but it is of course understood that the integrals extend to infinity in 
the appropriate direction. 

In addition to the differential identities (J3.10p ~ p.l4p for the functions 
Qt, Pt, and Rt , there are further identities in |241 125] that involve the four 
functions of the variable t defined by 

q{t) = Qt{t), p{t)=Pt{t), u{t) = Rt{t,t), v{t) = ^{u^-q^). (5.1) 

These four functions satisfy the closed differential system 

q = p — uq, p = tq + uq — 2vq, u = —q , v' = —pq. (5-2) 

The second order system ()3.2p gives formulas for -^tpi and ■^'ip2- We 
start by verifying the latter. 



From (j3.8p and (|3.9p we obtain 

q2 q 

_^2 = -2rlfi—tl)2 + {rlz + 2r2S2 - r2/i^)V'2 



/oo 
{x - t)F{z + C{x - t))Rt{x, t)dx. 
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Here we use (l2:20]l . (l3J3]l and ([O]) to obtain 

(x - t)Rt{x, t) = pQt{x) - qPt{x) = q^Rt{x, t) + q'Qt{x) - qQ'tix), 
where p = p{t). Thus by ([321) and dMD, 

Q'i Q 

_^2 = -2rjn—i)2 + ir\z + 2r2S2 - r|^^)V'2 

F(z + C(x - t)){q'^Rt{x, t) + q'Qt{x) - qQ't{x))dx 
= -2ri/i^^2 + {rlz + 2r2S2 - rf^^^^^ + Cr2g2(^2 _ F{z)) 

/oo 
F(z + C7(2;-t))Q;(x)(ix. 

We apply integration by parts to the remaining integral. The integrated 
term is Cr2q^F{z) and so we obtain (again using (|3.7p ) 

—■02 = -2r2/U— 0^2 + {rlz + 2r2S2 + Crlq^ - rl^?)ip2 

/oo 
F'{z + C{x-t))Qt{x)dx, 

d 

= -2rlfi—ij2 + (rlz + 2r2S2 + Crlq^ - r\ix^)^2 

- C^rlq—Tpi - CrlqTpi 

which is what is needed for -^ip2 according to (13. 2p . 

The proof for ■§^4'i is a bit more involved. We have by (|3.7p and p.9p 

q2 j-oo 

Q^^i = - F"{z + C{x-t))Qt{x)dx 

d 

= -2rln—i)i + {rlz + 2r2S2 - r2fi'^ - Crlt)il)i 



-Crl xF{z + C{x-t))Qt{x)dx. (5.3) 

Here we use p.l4p and (j3.13p to obtain 

xQt{x) = 2vQt{x) -pRt{x,t) + P;{x) - uPt{x) 

= -{u^- 2v)Qt{x) -{p- uq)Rt{x, t) + P/(x) - uQ't{x), 
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which by (|5.2p and the formula for v{t) in (j5.ip gives 

xQt{x) = -q^Qt{x) - qRt{x,t) + ^/(x) - 'uQj(x). 
Using this in (|5.3p we find using the definitions (j3.7p and (|3.8p . 

/oo 
F(z + C{x - t)){P[{x) - uQ[{x))dx. 

We integrate by parts on the remaining integral. The term —Cr2q'F(z) is 
cancelled by the integrated term Cr2{p — uq)F{z) because of (15.2p . Thus 

^2 Q 

-^i)i = -2r2/x— V'l + {rlz + 2r2S2 - 4li'^ - Cr^t - Crlq^)^^^ 

F'{z + C{x - t)){Pt{x) - uQt{x))dx. 
Now we use (I3.13P again and we simplify using (13.8P 

Q'i Q 

^1 = -2rl^i—^i + {rlz + 2r2S2 - 4l? - Cr^t - Crlq^)il^i 
OZ'^ oz 

/oo 
F'{z + C{x - t)){qRt{x,t) - Q't{x))dx 

= -2ri/x— V'l + {rlz + 2r2S2 - 4l^^ - Crlt - Crlq^)i>i 

+ Crl(i^2 - C\lq ("^^2 - F'{z)\ + C\l J" F'{z + C{x - t))Q[{x)dx. 

We integrate by parts, the integrated term cancels with C^r2qF'{z) and so 

_^i = _2r2/i— V^i + {rlz + 2r2S2 - r^fi^ - Crlt - Crlq^)^i 

- CVig— V2 + Crlq'i^2 - C^l / F"{z + C{x - t))Qt{x)dx. 

The final integral we can express in terms of -^ipi (see the first identity in 
(|5.3p ) and the result is that 

q2 q 

-^Vi = -2rliJ—i>i + {rlz + 2r2S2 - r^p? - Crlt - Crlq^)i)i 



- C\lq^i;2 + Crlq'ij2 + C^l^. 

oz OZ' 



26 



Rearranging terms we find 

- Crl(iil>2 + {-rlz - 2r2S2 + 4^i^ + Crlt + Crfg^)^,. 
Since C^''r2 — 1 = {r2/rif', this is 

2 

- CrfgVa + (-r-?z - 2^S2 + r^rf^^ ^ ^^2^ ^ Cr^g^)^^. 

^2 



Since t = ^-^ (2^ + 2ff - (^I + ''2)/^ J > we finally obtain 

^2 d d 

—^1 = cVfg— V'2 + 2r^/x— ■01 

- Crlq^2 + {-rlz + 2risi - rf^^ + Crjq^)'ipi 
as claimed in (13.211. D 
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